A SOLUTION TO THE PROBLEM OF THE SKIN EFFECT WITH A 
DISPLACEMENT CURRENT IN THE MAXWELL PLASMA BY THE METHOD 

OF EXPANSION IN EIGENFUNCTIONS 
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A problem of the skin effect in the Maxwell plasma is solved analytically by the method of 
expansion in eigenfunctions based on the Vlasov-Maxwell kinetic equation with a self-consistent 
electric field. Specular electron reflection from the boundary is used as a boundary condition. 
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I. INTRODUCTION 



The skin effect is caused by the electron gas response 
to an external variable electromagnetic field tangential 
to the surface [1]. This classical problem has been stud- 
ied by many authors (for example, see [13]). The present 
work develops an analytical method of solving boundary 
problems for systems of equations describing the behav- 
ior of electrons and an electric field in the half-space of 
weakly ionized plasma. This method is extremely con- 
venient, because it allows the sought-after distribution 
function to be derived in an explicit form. The method 
being developed is based on the idea of expansion of the 
solution in generalized singular eigenfunctions of the cor- 
responding characteristic system [2] obtained after vari- 
able separation. A solution to the characteristic system 
in the space of generalized functions [4] gives eigenfunc- 
tions with a continuous spectrum covering the entire pos- 
itive real semiaxis. The structure of the discrete spec- 
trum is elucidated by finding zeros of the dispersion func- 
tion, and eigenfunctions of this spectrum are determined. 
A general solution to the system of the Vlasov-Maxwell 
equations is constructed based on solutions for continu- 
ous and discrete spectra. The proof of the expansion in 
the eigenfunctions is reduced to a solution of the integral 
equation with the Cauchy kernels. The last is reduced 
to the Ricmann boundary problem in the theory of func- 
tions of complex variables. The solvability conditions 
and the Sokhotskii formulas allow all unknown expansion 
coefficients in the solution of the initial boundary prob- 
lem to be calculated. Let us assume that the Maxwell 
plasma occupies the half-space x > , where x is the 
coordinate orthogonal to the plasma boundary. Let the 
external electric field has only one y-component. Then 
the self-consistent electric field inside the plasma will also 
have only one y-component E(x)e~ lujt . We now con- 
sider the kinetic equation for the electron distribution 



function: 

| + v,| + e^)e-^ = ,(/o-/(M,v)). (1) 

where v is the frequency of electron collisions with ions, 
e is the electron charge, and fo{v) is the Maxwell equi- 
librium distribution function: 



3/2 



/o(v) = n(£) exp(-/3V), (3 



2k B T 



Here k is the Boltzmann constant, T is the plasma tem- 
perature, v is the electron velocity, m is the electron 
mass, and n is the electron concentration. 

Let us assume that the field strength is such that the 
linear approximation is applicable. Then the distribution 
function can be represented in the form 

/ = /o (1 + C y exp(-iu}t)h(x, n)) , 

where C = \f]3v is the dimensionless velocity of electron 
and /x = C x . We now introduce dimensionless quantities 
t\ = vt, x\ — v^f]3x , and 



e(xi) 



V2e 
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E( Xl ). 



Then we will write again x instead of x± .In the new 
variables, kinetic equation (1) and the field equation with 
allowance for the displacement current are written as fol- 
lows: 



9h , , , , , 
^dx + Z ° = 



(2) 



e"(x) + Q 2 e(x) 



cxp(— ji ) h(x, fjf) d/jf , (3) 



c 



z = l 



ILOT, 



2 



where I is the free path of the electron, S = 



2ttuj<to 



is the classical depth of the skin layer, oq = , cto 

mv 

2l 2 

is the electric conductance, a = —r- , a is the anomaly 

o z 

parameter. 

Let us formulate conditions for the distribution func- 
tion and field on the plasma boundary: 

h(0,n) = h(0, -n), 0<m<+oo, e(0) = e s . (4) 

We search for a distribution function and field that 
decay with increasing distance from the surface: 

h(+oo, fi) = 0, — oo < ii < +oo, e(oo) = 0. (5) 

Without loss of generality, we further set e s = 1 . 

II. EIGENFUNCTIONS AND EIGENVALUES 



Separation of variables (see [2] ) 

x 

h v (x, jj,) = exp(-z -)$(?y, /u), 



e v (x) = exp(-z -)E(ri), 

where r\ is a complex spectral parameter, reduces system 
of equations (2) and (3) to the characteristic system 



where 



[^ + QV]%) = ->(n), 
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i(v) = J e^ 2 $(??,^)a>. 



(6) 



(7) 



From Eqs. (6) and (7) we find the eigcnfunctions of the 
continuous spectrum in the class of generalized functions 
[3]: 

M ) = "n^P— + \(n)6( V M ), (8) 

y/TT T] - H 
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(9) 



Taking into account the decrease of the distribution 
function and electric field far from the boundary, the pos- 
itive real semiaxis < x < +oo is taken to mean the con- 
tinuous spectrum of the boundary problem. The eigen- 
functions of the continuous spectrum h n (x,iJ,) and e n (x) 



are decreasing functions of the variable x for 3?zo > . 
The eigenf unctions in equalities (8) and (9) have been 
normalized by the condition 



oc 

/ 



and the dispersion function 
X(z) = 1 + 



1 + 



z* + 



az 

7^ 



ji — z 



has been introduced. 
Let us designate b 



Q 2 



and express the dispersion 



function of the problem in terms of the dispersion func- 
tion of the Van Kampen plasma Xo(z) : 



X(z) = l + (b- a)z 2 + az 2 X {z), 



where 



Ao(*) 



fie M d[i 

[1 — z 



For the dispersion function in the vicinity of the point 
at infinity, the asymptotic expansion 



X(z) = (b- af + (!--) 



3a 
4? 



15a 



oo, 



is fulfilled. 

We now elucidate the structure of the discrete spec- 
trum by the method developed in [2, 3]. By definition, 
this spectrum consists of zeros of the dispersion function 
laying outside of the cut (— oo, oo) . 

Let N be the number of zeros. Since the dispersion 
function has a double pole at the point z — oo , the 
number of its zeros is 



^ = 2+^[argA(z)] 7e , 



(10) 



where j e is a contour passing clockwise over the cut 
(—00,00) at distance e and having no zeros inside. 
Taking the limit in Eq. (10) when e — > , we obtain 



N = 2+ 



1 

2^ 



arg 



A + (r) 
A-(r) 



2+- 



( — 00,00) 



arg 



A+(r) 



A-(r) 



(0,00) 



Here A ± (t) = A(/x) ± iira^e^^ 2 are the maximum 
and minimum values of the function X(z) in the cut. 

Let us consider the region D + (we designate by D~ 
its external boundary) in the a plane whose boundary 
is set by the parametric equations 

dD+ = {a = ai + ia 2 : 3?A+(^) = 0, 



3 



3?A+(/j) = 0, -oo < [i < oo}. 

By analogy with [2], we can demonstrate that 1) if 
a e D+,N = 4 and 2) if a e D~,N = 2. The mode 
with a G is not considered here, since it has already 
been studied in detail in [3]. 

Let us write down (discrete) eigenfunctions corre- 
sponding to the obtained discrete spectrum {±77^ : 
Afe) =0, ft = 0,1}: 

* = ' 



k=0 



E{rt k ) = — i:: 



E 

fc=0 



(ft = 0,1). 



We note that in the last formulas, ft = when a e D 
and ft = 0, 1 when a G L> + . 



I 1 

-= V Afcjfi! exp(-^)+ 



k=0 



1 1 

+ V r/ 2 cxp(-7 ? 2 )A(7 7 ) ( i?7 = . (14) 

^k~o az ° 



where 



<P{p) = A=y^A k rf k cM-vl){— —)■ 

Let us transform Eq. (14) setting A(—rj) = —A(rf) , 
that is, expanding the coefficient A(rf) to the entire real 
axis as an odd one. Considering that $(—77,—^) = 
$(77,/^) , we reduce Eq. (13) to the form 

00 

<p(/j,)+ J A{rj)${ri,ii)dri = 0, -00 < \i < 00, 



III. ANALYTICAL PROBLEM SOLUTION 

Let us represent the general solution of system (2)- 
(5) in the form of expansion in eigenfunctions of the 
discrete and continuous spectra, automatically satisfying 
the boundary conditions at infinity: 



or after substitution of the eigenfunctions into this equa- 
tion, 



a f ij 3 A(ri) exp(— Ty 2 ) drj 



+ \(ji)A(n)+ 



^>-£££?W-rf-f) + 



+ 



OO 

y"exp(-^)^)$(^,/i)d»7, (11) 



+a(p(/j.) — 0, —00 < /j < 00. 
Let us introduce the auxiliary function 



N(z 



1 f T] 3 cxp(-ry 2 )A(?7) 



77 — z 



dij, 



(15) 



e(x) = — > A^exp 



Vk 



whose boundary values, according to the Sokhotskii for- 
mulas, obey the equality 

N+{p) - N-(n) = 2V¥^ 3 cxp(-7i 3 )A(7i) = 



OO 

+ ^ jexp(- t, 2 -^) 7^(77)^. (12) 


Here A k (ft = 0, 1) arc unknown coefficients of the 
discrete spectrum with A\ = for a e L>~ , ^.(77) is 
unknown function called the coefficient of the continuous 
spectrum, 3i(z n /rjk) > (ft = 0, 1) , and 3f?z = 1 . 

Substituting expansions (11) and (12) into the bound- 
ary conditions, we obtain the following integral equa- 
tions: 

OC OO 

aip{p)+J A^)^,^- J A( V )$(r 1 ,-n)dri = 0, (13) 



= ^ [A +0,)_ A" (/*)]. 
a 

With the help of boundary values of the auxiliary func- 
tion N(z) and the dispersion function, we reduce the in- 
tegral equation with Cauchys kernel (15) to the Riemann 
boundary problem 

A + (/i)[J\r+(/i) + ip(n)] - X-{im)[N-(jm) + ^)], 

whose general solution has the form 

! 1 

N(z) = -—^2 A kV l cxp(-77 2 )x 

* k=0 
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Vk 



Chz 

wr 



(16) 



Eliminating the first-order poles at points rj k , we ob- 
tain 

Ci =-^=A kV ieM-vl)X(r 1 k) (fc = 0,l). 

Substituting general solution (16) into the Sokhotskii 
formula, we obtain the coefficient for the continuous spec- 
trum: 



rf exp(-77 2 ),4fa) 



2^i 



\+(rj) A" fa) 



We now return to Eq. (14) and write it in the form 



A' fan) A 



oo 

i , j_ r 

'(77O + 2m J 



A+fa) A" fa) 



dry 



azoCi 



(17) 



After integration of Eq. (17) by the methods of contour 
integration, we transform the last equation and calculate 
first the constant C\ : 



oo oo 

Cl = azoJ(ay J{a) = 2^ / Mir) = nj 



dr 



and then constants A k with the help of Eq. (14): Ak = 

V^ C M4) (fc = ,i). 

az J(a)r] 2 k X'{r] k ) 

To calculate the impedance, we consider the electric 
field derivative 

e'(0) = azlCv 



Taking into account the evenness of the integrand, we 
obtain 



1 

2^i 



1 



1 



A+fa) A" fa) 



1 



1 



1 



dr] 

T] 2 rj \'(r) ) r) \'(r]o)' 



Now it is clear that the derivative of the electric field 



is e'(0) = ^fC, 



2J(a) 



and the expression for the 



surface impedance is 

„ 8iriuil 



c 2 z 



dr 

it J A(zt) 





(19) 



Let us express all constants in Eq. (19) in terms of 

u v Anne?. . 
7 = — and e = — , where lo p = is the plasma 



UJ p 

frequency, b 
l 



p 
7 2 



m 



(e - i'y) 



V r , a = — i 



7 



2 c ' 



(e — «7) 3 



v 2 , and 



v c = 




Fig. 1. Modulus of the impedance. 
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+ 



%A'(r? ) r/iA 



OO 

1 — [ 

Vfai) 2wi J 



1 



1 



A+fa) A" fa) 



dr] 



To integrate this expression, we take advantage of the 
representation 



A( 



- = — / 

(z) 2m J 



1 



A+fa) A" fa) 



dr] 

r] — z 



E (n? - 



2rjk 



^fa 2 -z2)A'fa fc )' 
From equality (18) for z = , we obtain 



£ 

k=0 



+ 



r] k X'(r] k ) 2iri 



OC 

- / 

Im J 



(18) 



A+fa) A- fa) 



dr] 



Fig. 2. Argument of the impedance. 



Let us now represent dispersion function (10) in the 
form 

X(z) = 1 + - l2y2 \, z 2 + i . n(z) = 



(e — ij) 2 (e — ij) 2 



5 



l Z ol 



tude of the normal skin effect and Zq is the dimen- 
sionless impedance. The behavior of the dimensionless 
impedance modulus is shown in Figs. 1 and 3, and the 
behavior of its argument is illustrated by Figs. 2 and 4 for 
e = 1(T 3 and v c = 1CT 3 . The plots in Figs. 3 and 4 are 
drawn near the resonance, that is, when the parameter 
7 passes through the value 7 = 1 for u) = lo p . 



Fig. 3. Modulus of the impedance. 



An analysis of plots drawn in Figs. 1-4 demon- 
strates that near the plasma resonance, the modulus of 
the impedance has a sharp maximum which is not ob- 
served in the low-frequency limit or in the theory of 
normal skin effect, and the argument of the impedance 
changes abruptly near the resonance. 



Fig. 4. Argument of the impedance. 



(e — ij) 3 



[(e — ij) 3 + (e — ij)~f 2 v 2 z 2 + i^v 2 c z 2 + ijv 2 p(z)] , 

-.3 oo £xp^ — M 1 

where p(z) = = f d/i. After substitution 

of the expression obtained into formula (19) for the sur- 
face impedance, we have 

Z = — = x 

& Zq 



(e — ij) 3 dt 



n J — ij) 3 + (e — i"f)j 2 v^t 2 + i^v 2 p{t) 
o 



IV. CONCLUSIONS 

The expression for the impedance can be represented 
as Z = RZq , where R = 2itlu5c~ 2 is the magni- 
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